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NONEXISTENCE OF STEADY SOLUTIONS FOR ROTATIONAL SLENDER 
FIBRE SPINNING WITH SURFACE TENSION 

THOMAS GOTZ AND AXEL KLAR 


Abstract. Reduced one-dimensional equations for the stationary, isothermal rotational spin¬ 
ning process of slender fibers are considered for the case of large Reynolds (5 = 3/Re <C 1) and 
small Rossby numbers (e <C 1). Surface tension is included in the model using the parameter 
K, = yjlxI (2We) related to the inverse Weber number. The inviscid case 5 = 0 is discussed 
as a reference case. For the viscous case 5 > 0 numerical simulations indicate, that for a 
certain parameter range, no physically relevant solution may exist. Transferring properties 
of the inviscid limit to the viscous case, analytical bounds for the initial viscous stress of the 
fiber are obtained. A good agreement with the numerical results is found. These bounds give 
strong evidence, that for 5 > (l — -h no physical relevant stationary solution can 

exist. 

Keywords: Rotational Fiber Spinning, Viscous Fibers, Surface Tension, Boundary Value Prob¬ 
lem, Existence of Solutions 


1. Introduction 


Being an important technology, fiber spinning has recently gained considerable attention in the 
mathematical literature, see El El n El 0 El 0 0 uni HU [H ng . Describing a slender fluid jet, 
the model equations covering spinning processes are settled on balance laws for mass, momentum 
and energy. Usually, due to slenderness of the fiber, cross-sectional averaging is performed to 
obtain spatially one-dimensional models [lEIlilEli- 


In the present work we consider the isothermal, rotational spinning of slender fibers. As an 
example of industrial relevance, one might have the production of glass fibers in mind. In this 
situation energy balance holds and Coriolis and centrifugal forces as well as surface tension effects 
have to be included, see [210E3]. Neglecting gravity the reduced model equations read in 


as 


(l.la) 

dfA + ds{Au) = 0 , 

( 1 -lb) 

dtv + udgV = 5— 

( 1 . 1 c) 

5(7 -1- u5s7 = V , 

(l.ld) 

l|a«7ll = 1 • 


ds{Ads^dsu) ds{\fAds^) 


A 


A 


+ -V + 


The independent variables are the arc-length s € [0, L] and the time t G [0, oo). By 7 ( 5 , t) G 
and v{s,t) G we denote the position and the velocity of a point located on the fiber’s cen¬ 
terline at time t. Additionally, u{s,t) G R describes the tangential velocity of a fluid particle 
moving along the centerline. By A = A(s, t) we denote the cross-sectional area of the fiber. 
The parameter S = 3/Re is related to the Reynolds number Re and e is the Rossby number, 
i.e. the inverse of the scaled rotation frequency. The surface tension parameter k = yT/(2We) 
is related to the inverse Weber number [Tg. Equation (ll.Ial) is the continuity equation. The 
momentum equation (ll.lbl) includes surface tension and viscous, Coriolis and centrifugal forces. 
Equation (ll.lcl) describes the fact that the total derivative of 7 is given by v and finally. Equa¬ 
tion (jl.Idjl is the condition for the arc length parametrization of the fiber which is defined by a 
normalized tangent. 
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We mention that equations (HID have been derived from the full set of three-dimensional Navier- 
Stokes equations by an asymptotic analysis in the above mentioned papers. It is understood 
in many papers that in a boundary layer close to the orifice the above slender models do not 
describe the jet well. In this small region near the orifice the flow is fully three-dimensional and 
other effects might become important. In the present work we identify parameter regions where 
the slender model breaks down and can not be used any more to describe fiber spinning. For 
further investigation on the breakup of liquid jets from a rotating orifice we refer to [l4l [T5] . 

In application relevant cases one is typically interested in stationary, i.e. time independent solu¬ 
tions for small values of the parameters. The purpose of our present investigations is to identify 
those parameter regimes for 6, e and k where physically relevant stationary solutions exist for 
equations (HD . In a previous work [16] we analyzed the case without surface tension, i.e. k = 0. 
Applying similar techniques, we are able to extend these results to the model including surface 
tension. The results will show, that surface tension does not help to overcome the problem of 
non-existence of physically relevant solutions — a hope that was expressed in the outlook of our 
earlier work |16] . 

Finally, we mention that in m a one-dimensional rod model allowing for stretching, bending 
and twisting has been investigated asymptotically and numerically, which allows for simulations 
in the whole range of parameters. Moreover, a more detailed analysis of the parameter regions 
where the string model breaks down is given there for the case without surface tension. 

This work is organized as follows: In Section [21 we rewrite the stationary version of the above 
system p.ll) in a more tractable form. Appropriate boundary conditions are specified. Section|3] 
is devoted to the inviscid case <5 = 0. The main result is contained in Section |4l where we 
characterize the stationary solutions of the viscous problem <5 > 0 and show the non-existence 
of physically reasonable stationary solution for some range of the parameters <5, e and k. In 
Section [5] we present numerical simulations based on a solution of the boundary value problem 
in Eulerian coordinates. The obtained results confirm the analysis. 


2. The stationary case 


In Eqn. (Il.lal) . the mass flux is given by Au. In the stationary case, the mass flux is constant 
and after appropriate scaling one gets Au = 1. Physical meaningful solutions are characterized 
by strictly positive it, hence A = 1/u. Furthermore, (ll.lcl) reads as v = udg^ and inserting this 
into (ll.Ibll yields 


(2.1a) 

(2.1b) 


ds {u dsj) = Sds 


^dsud^y^ 


115.711 = 1 • 



2 a J. 11 

-k -<9.7^ + —-7 , 
e u 


To eliminate the algebraic constraint (I2.1bll . we project dgj onto the unit circle and parameterize 
(?s7 := T = (cos Q!, sin q;) by the angle a = a(s). Furthermore, we use polar coordinates 7 = 
r(cos</5, sini^) for the centerline and introduce jS = a — tp. Abbreviating derivatives with respect 
to s by a prime we transform (|2.1ap via q = u — 5u'/u — nj^/u into a first order system which 
we supply with boundary conditions at s = 0 (nozzle) and s = L (end point of the fiber): 


(2.2a) 

(2.2b) 
(2.2c) 
(2.2d) 


Su' = u { u -;= — q 

y/u 


e^uq' = r cos 13 , 


r' = cos/5 , 

P' = --- 

eq 


e^uq 


+ 1 


sin/3 


u(0) = 1 , 


9k 

qiL)=u{L)-^= 

y/u{L) 

r(0) = 1 , 


m = 0 . 


Remark 2.1. The above boundary conditions for u, r and /3 given at s = 0 express the fact 
that the fiber leaves the nozzle tangentially with a fixed velocity. Due to the boundary condition 
for q at the end of the fiber, i.e. at s = L, the viscous stresses at the end of the fiber balance 
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the surface tension. Note, that in the inviscid case d = 0 without surface tension, i.e. k = 0, we 
have due to positive u the equation u(s) = q{s) for all s. 

Lemma 2.2. In (j2.2p it holds that q'{0) > 0. 

Remark 2.3. The observable q can be interpreted as the internal energy of the fluid being 
the sum of the kinetic energy Au^ = u and the viscous stress energy —SAu' = —Su'/u and the 
potential energy —ny/A = —Kj^/u due to surface tension. 

Remark 2.4. The polar angle (p of the centerline is determined via 

(2.3) <p' = ^ , (^(0) = 0 

r 

leading to <^'( 0 ) = 0 . 

Remark 2.5 (Equations in Lagrangian coordinates). For describing the parameter regimes 
for 6, e and k where physically relevant solutions of ( 12 . 21 ) exist, it is convenient to introduce 
Lagrangian coordinates. The passage from Eulerian to Lagrangian coordinates is settled on the 
transformation 

= u{s(t)) > 0 , s( 0 ) = 0 

mapping the Euler position s of a material point to the Lagrangian parameter t. Starting from 
the nozzle it takes a fluid particle t time units to approach position 7 ( 5 ). It takes a particle T 
time units to travel from the nozzle to the endpoint s = L. 


After some elementary manipulations we obtain the Lagrangian system 


(2.4a) SiiL = u\[ul- 


y/uL 


dL 


(2.4b) e^qL = tl cos (II , 

(2.4c) t-l = UL cos Pl, 

(2.4d) p^ = -^- 

eqL 


s^qL 


UL 


sin/Si 


tl 


ul{0) = 1 , 

qL{T) = UL{T)- 
r-L(0) = l , 

/3l(0) = 0. 


2k 




Let UL(t) = u{s(t)) be the velocity in the Lagrangian framework and mutatis mutandis for the 
other variables. By a dot, e.g. 

duL{t) du(s(t)) ds(t) . 

^ dt ^ Js ^ ^ ’ 

we denote the derivative with respect to the Lagrangian parameter. 


3. Inviscid limit (5 = 0 


As outlined in the Introduction, we are interested in cases where (5, e and k are small. In this 
Section we shall exploit this fact in a formal manner. As small e corresponds to a high rotation 
frequency one expects u = 0(l/e) at least in the fiber’s interior. Due to the definition of q one 
also expects q = 0(l/e). Therefore it is convenient to introduce the re-scaled variables 




V = su , 

w = eq . 

In terms of 

v, w, r 

and P the model equations ( 12 . 21 ) become 

(3.1a) 

dev' 


v{0) 

(3.1b) 

vw' 

= r cos p , 

v{L) 

(3.1c) 

r 

= cos P , 

r( 0 ) 

(3.1d) 

P' 

2 / sin ,5 

w V vw J r ’ 

m 


e , 


w(L) + 


A 

y/v{L) ’ 


1 , 


0 , 
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Figure 3.1. Phasepor- 
trait of the system in 
the r/3-plane. 


Figure 3.2. Velocity 
V = eu oi an inviscid fiber 
with and without surface 
tension. 


where A = In the —formal— inviscid limit d = 0 the system (|3.1I) together with the 

assumed positivity of v yields 

Lemma 3.1. In the inviscid limit S = 0 of the stationary system (EH) it holds that 

X K 

w = v - — i.e. q = u - 

y/v y/u 


In the inviscid limit <5 = 0, we can eliminate w in the above system EH and obtain the reduced 
system 

-1 


(3.2a) 

(3.2b) 

(3.2c) 


A 

‘iy/v 


r cos (3 , 


r' = cos /3 
/?' = - 


V - j= \v‘^ - \y/v 

y/v ^ ^ 


+ 1 


sin/3 


u(0) = e , 
r(0) = 1 , 

m = 0 . 


This system can be solved numerically. The Figures ITTl and lT^ show simulations for e = 0.16 and 
K = 0.5 leading to A = 0.032 compared to the results without surface tension, i.e. /c = A = 0. The 
numerical results show, that including surface tension to the inviscid model alters the behaviour 
of the fiber only marginally. 


Lemma 3.2. For 0 < e < — it holds that 


(3.3) 


1 + k/2 

(I -L n - I 

£^1 + !)' 


Proof. Considering the system EH, the expression for v” reads as 


V = 


r' cos /3 — r/3' sin /3 , ^ ^ 

-r- rv cosp- 


2 v^ 




and using the boundary conditions at s = 0, we obtain 


v"(0) = 


£3(l + f)^ 


□ 
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Lemma 13.21 shows, that in the inviscid, fast rotating limit the velocity in Eulerian coordinates is 
locally at the nozzle a concave function. 

Lemma 3.3. For the second derivative of the Lagrangian velocity ul it holds that 
(3.4) ilLiO) > 0 . 


Proof. Using (12.41) with 5 = 0, we get 


e^iiL = tlcos/^l 1 + 


3/2 


2u^ 

Differentiating again and inserting the conditions at s = 0 yields 

1 + f) 


wl(0) = 


> 0 . 


□ 


Lemma |3.3I shows, that in the Lagrangian framework of the inviscid limit 5 = 0, the velocity is 
a convex function at least close to the nozzle. Note, that ii > 0 implies that the acceleration of 
a fluid particle increases with its flight time. This agrees with the fact, that the acting forces, 
i.e. the centrifugal and Coriolis forces, also increase the further outward the fluid particle travels. 


4. The viscous case 


After the investigation of the inviscid case 5 = 0, we turn our attention to the boundary value 
problem (BVP) (12.2|) governing the stationary viscous case 


(4.1a) 


5m' = u 




(4.1b) e^uq = r cos /3 , 


(4.1c) 

(4.1d) 


r' = cos /3 , 



m(0) = 1 , 


q[L) = u{L) 
r(0) = 1 , 

m = 0 . 


2k 

\/u{L) ’ 


To gain some insight into the behaviour of the solutions for small values of the parameters 5, e 
and K, we carry out numerical simulations. The results are obtained with Matlab® 2008a using 
the routine bvp5c, a finite difference code that implements the four-stage Lobatto Ilia formula. 
This is a collocation formula and the collocation polynomial provides a -continuous solution 
that is fifth-order accurate uniformly on the considered interval. The formula is implemented 
as an implicit Runge-Kutta formula and solves the algebraic equations directly. As an initial 
guess for the solution we use the inviscid case 5 = 0. 


5 

e = 0.25, K = 0.1 
5 (0) m"(0) 

m(0) 

K 

£ = 0.1,5 = 0.01 
5(0) m"(0) 

m(0) 

0.1 

0.1252 

-4.8225 

55.221 

0.3 

0.1535 

-208.43 

2797.9 

0.125 

0.0239 

-3.7310 

45.307 

0.4 

0.0653 

-194.37 

2681.0 

0.13 

0.0092 

-3.5535 

43.348 

0.45 

0.0220 

-168.59 

2617.5 

0.133 

0.0021 

-2.2389 

43.288 

0.475 

0.0013 

-149.47 

2592.1 

0.135 

no 

convergence 

0.48 

no 

convergence 


Table 4.1. Results of the numerical simulations using Matlab®. 


The numerical results, some of them are given in Table 14.11 as well as the analytical results 
for the inviscid case (see Lemmata 13.1113.21 and [T3l) show, that a solution of the stationary 
system m, if it exists at all, has the following properties: 
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Definition 1 (Physically relevant stationary solutions). We call a solution (u,q^r,P) of the 
boundary value problem (ED a physically relevant stationary solution, if it has the following 
properties 

(PI) go := 9(0) G (0,1 - «] , 

(P2) u"(0) < 0 , 

(P3) ul(0) > 0 . 

Remark 4.1. In the sense of this definition, all the numerical solutions we were able to compute 
are physically relevant; the properties (iPip. (IP^ and (|P3ll are satiesfied. Furthermore, the 
stationary solution of the inviscid case, i.e. the solution of (13.2p together with the according 
solutions for u and q is physically relevant. 

Let us discuss the properties (IpD)—(IP 3D in detail. 

Lemma 4.2. Let {u,q,r,fi) be a solution of the BVP (14.ip . Then, property (IPlIl is equivalent 
to u'{0) > 0 and /3'(0) < 0. 


Proof. Eqn. (I4.1al) reads as 


ufO) 



1- K-qo 


S 


and hence 


m '( 0 ) > 0 <;=^ qo ^ ^ — K . 

Due to (|4.1dp 



holds and /3'(0) < 0 if and only if go > 0. □ 


Remark 4.3. The conditions u'(0) > 0 and ,8^(0) < 0 allow for an easy interpretation in physical 
terms: 

The condition u'(0) > 0 assures, that the fiber is not accelerated into the nozzle—a fact that 
one should expect from a physically relevant solution. 

Due to R,emark [T4l we observe that the condition /3'(0) < 0 is equivalent to a'(0) = ,8'(0)—V5'(0) = 
,8'(0) < 0. Hence the fiber initially bends in the sense of rotation of the drum. The opposite 
condition q;'(0 ) > 0 would imply, that the fiber initially bends against the direction of rotation 
of the drum. This is obviously counterintuitive and physically not reasonable. 

Remark 4.4. The property (IPlIl also implies, that physically relevant solutions require 0 < k < 

1 . 


Lemma 4.5. Let {u,q,r,/3) be a stationary solution of the BVP (14.11) . Then, the property (IPID 
implies 

(5V'(0)< (2-0 (!-«)-A. 

In other words: For all 6 > 0 and -r—--; (IPlI) (IP2I1 . 

( 2 -§)(!-«) 


Proof. Via differentiation we deduce from (14.ip the equation 


S^u" = u' { u — — q ) +u I u' + 


‘^\/u 


- 9 ' I = T ( 2m - 


2 -v/m 


q\\u 




rcos/3 


As a consequence we obtain at the nozzle, i.e. for s = 0 

(4.2) d’^u'fO) = (2 - - go01 - K - go) - ^ , 
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and the property m yields 

<5V(O)<(2-0 

□ 


7 -7—-As this statement 

( 2 -§)(!-«) 

holds for fixed 6 for all “small” parameters e, we can deduce (IP2I) rigorously from m in the 
case of small e and small k. In particular, for fixed 5, the property (IP2I) does not result in a 
restiction on the values of £ and k as £ and k tend to 0. 


This results shows, that u"{Q) < 0 for all 5 and £ with 2£^ < 


The property (IP3I) of a physically relevant solution is required due to numerical evidence and 
analogy to the inviscid case,. See also the arguments at the end of section[31 It will lead us now 
to bounds for and finally to conditions on k, e and 5. 

Theorem 4.6. Let (u,q,r,l3) be a physically relevant stationary solution of the BVP (14.11) . 
Then, the following estimate for the initial value qo holds: 


where 


max (O,^) <qo< min (1, qo) 


3-f+ 

qo = - 2 


Proof. Due to (j4.2l) we have 

S’^u'fO) = (^2 - - go) (1 - « - 9o) - ^ = M • 

The property (ira requires /x < 0 and hence we obtain 

3-¥-+ 3-f+ v/(l + f)'+4fe-2 

90 =-^- < 90 < -^-> 1 ■ 

Analogously, due to 1x^(0) = xx"(0)u(0)^ + xx'(0)^it(0) and using (14.21) again, we get 


S^ul{0) 

Now, (IP3I) . i.e. A > 0 leads to 


(1 - K - 9o) ( 3 - -K - 2qo 



5-^ + ^(l + f)V8&-2 5-f-^(l + f)V8fe-2 

90 > - ^ - >1 or 90 < 90 = -- 

Summarizing the individual estimates shows the assertion. 

Corollary 4.7 (Non-existence of physically relevant solutions). If 

(4-3) p { k ) = 3(^i-Y + Yj<^ 

no physically relevant stationary solution exists. 


□ 


Proof. If 9 o = 0, physically relevant stationary solutions cannot exist. 


□ 


Remark 4.8. The results reported in Table l4^ confirm these findings The boundary between 
numerical convergence and not-convergence is close to the analytical non-existence result (14.31) . 
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(5 

5 

n 


p{k) 


0.133 

0.25 

0.1 

2.128 

2.565 

convergence 

0.135 

0.25 

0.1 

2.16 

2.565 

no convergence 

0.01 

0.1 

0.475 

1 

1.201 

convergence 

0.01 

0.1 

0.48 

1 

1.186 

no convergenee 


Table 4.2. Numerical results for different values of the parameters 


Remark 4.9. The result of the corollary can be partly understood, if the parameters are rein¬ 
terpreted in their physical context. We have, that <5 = 3/Re is proportional to the kinematic 
viscosity v. Furthermore e, being the Rossby-number, is proportional to the inverse of the ro¬ 
tation frequency w of the spinning drum and k is proportional to the surface tension a. With 
this interpretation in mind, Corollary 14.71 states, that with some scaling constant c > 0 and with 
p(u) = p{k) , we have that 

Q 

• ioT v < V := —p(a) physically relevant stationary solutions may exist and 

_ ojf 

• for V > V physically relevant stationary solutions cannot exist. 

Hence, for a given rotation frequency u and surface tension cr, fibers with a viscosity v > v 
cannot be spun. The larger the surface tension, the smaller the limiting viscosity V has to be. 
Moreover, there is a limit surface tension ct related to the value k = 1. In case a > a there is no 
solution, independently of the values of e and S. 

5. Numerical results 

The Figures [3T] and [321 on page [4] show simulations for e = 0.16 and k = 0.5 compared to the 
results without surface tension, i.e. k = 0. These simulations indicate that including surface 
tension alters the inviscid case just marginally. 

Figure [5T] shows the parameter values, for which the equations (14.11) were solved successfully. 
Table |4T] reports some of these results. The solid line corresponds to the theoretical bound (14.31) . 
In Corollary 14.71 it was shown that only in the region below the solid line physically relevant 
stationary solutions may exist. Our numerical simulations confirm this finding. 

Figure 112] illustrates the behavior of qo for fixed surface tension k = 0.1. Simulations are shown 
for e = 0.1 and e = 0.24 as well as for different values of 5. The symbols indicate the numerical 
results, whereas the lines show the bounds for qq due to Theorem 14.61 It is clearly visible that all 
our numerical results are within the limits predicted due to the assumption of physically relevant 
solutions. 


6 . Conclusions 

There is strong numerical evidence that the regime where physically relevant stationary solutions 
exist is characterized by the inequality 

( 6 . 1 ) 6 < e^p{K) . 

Naturally, one has to interpret a coupling of S, k and e like dSIl) in terms of the kinematic 
viscosity u of the fluid, the rotational frequency w of the drum and the surface tension a. In 
these variables the statement above is translated into the following: for v > v := —p(a) the 
fibre can not be spun. 

A detailed analysis—both theoretically and numerically— of the behavior of transient solutions 
to (HID may reveal, whether physically relevant transient solutions show up in the parameter 
regime, where stationary solutions cease to exist. 
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Figure 5.1. The param¬ 
eters (5/e^ and k for the 
successful simulations in¬ 
dicated by dots. The solid 
lines shows the theoretical 
bound (14.31) . 



Figure 5.2. Plot of qo 
vs. 6 for At = 0.1 and e = 
0.1, 0.25. 
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